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Abstract Dynes ansatz [T] is a broadly utilized procedure to extract the su- 
perconducting energy gap from the tunnelling differential conductance, a quan- 
tity that is measured by scanning tunnelling microscopy (STM) techniques. In 
this letter we investigate the limit of applicability of this ansatz in nano-scale 
superconductors and propose a generalization that permits to study thermal 
and quantum fluctuations in STM experiments. 

Keywords Superconductivity • STM • Dynes fitting • Thermal fluctuations 
PACS 74.20.Fg • 75.10.Jm 



1 Introduction 

Scanning tunneling microscopy techniques (STM) have revolutionized |2J the 
field of experimental superconductivity for both conventional and high tem- 
perature superconductors. A typical observable in these experiments is the 
normalized tunnelling differential conductance, 
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where N s is the local density of states, T is temperature, V is the gate voltage. 
In order to make contact with properties of the superconductor, such as the 
superconducting energy gap A, is required to fit the experimental data with 
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a theoretical expression of dl/dV that includes A as a fitting parameter. In 
principle this is not an easy task as for many materials the theoretical local 
density of states N s (uj, A) is not exactly known or its calculation involves the 
use of the rather cumbersome Eliashberg theory of superconductivity. However 
Dynes noted [T] that for some metals the energy gap A = An obtained by a 
simple fitting function, 

N D (.)=Im( . U + iF ) (2) 

is in excellent agreement with the one obtained by other experimental tech- 
niques. ro(T) in ([2]) is other fitting parameter that describes different sources 
of decoherencc. In light of this success this procedure is now employed to obtain 
the energy gap in STM experiments involving a broad range of conventional 
and high T c superconductors. In case of a material with a non s-wave order 
parameter symmetry the Dynes ansatz is modified to account for the expected 
momentum dependence of Ap. 

In this letter we investigate the limit of applicability of Q in nano-scale super- 
conductors where deviations from mean-field predictions are well documented 
(3j|4] . These were recently used in [5] to analyze STM data from Pb supercon- 
ducting nano particles. We focus on sizes L < lOnm for which deviations from 
mean-field predictions start to be relevant. First we introduce a model capable 
to describing thermal fluctuations and provide theoretical expressions for the 
gap and dl/dV. Then we fit this dl/dV by Dynes ansatz pi, fit and compare 
the resulting gap with (^.Finally we propose a generalization of Dynes ansatz 
that broadens the information obtained from STM experiments. 



2 Theoretical description of thermal fluctuations in zero 
dimensional superconductors 

In order to test Dynes procedure in the nanoscale we study the reduced BCS 
Hamiltonian by path integral techniques in the static path approximation 
(SPA) [6. which provide a good description of thermal fluctuations in the 
zero dimensional limit corresponding to grain sizes much smaller than the su- 
perconducting coherence length. Here we only present the main results and 
refer to 015] for details. The partition function Z is given by, 

Z = [d\A\ \A\ e-^-^CW) (3) 
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where £ = y e 2 + \ A\ 2 , /3 = 1/T and g (e) = J2 a S {e — e a ) is the one particle 
density of states, Zq is the parition function for non-interacting electrons, A 
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is the dimensionless coupling constant, Ed is the Debye energy and <5 is the 
mean level spacing. From the explicit knowledge of Z , the normalized density 
of states (DOS) N s (uj) is, 
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where i"s(T) is a phenomenological parameter that describes all decoherence 
mechanisms except thermal fluctuations. The superconducting gap A in this 
formalism is given by, 
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We note that A leads to the bulk gap A 
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for L — > oo however 



sinhl/A 

finite size effects in A might differ slightly from those in the spectral gap [3j. 



3 Justification of Dynes ansatz in the nanoscale region 

We now test Dynes ansatz by comparing Ajj(T) ^ with A(T) ^ where the 
former is obtained by fitting the theoretical dl/dV (obtained from (|5|) with 
the Dynes dl/dV coming from Dynes ansatz is presumably applicable pro- 
vided that the low energy excitations are well defined quasiparticles - Fermi 
liquid theory holds - and a gap do exist in the spectrum such that the quasipar- 
ticle dispersion relation is E = \/e 2 + A 2 even if A is not the BCS prediction. 
We show in Fig.JT^ that Dynes expression fits very well the theoretical dl/dV 
and in Fig. [TJd that the energy gap Ad(T), obtained from |2]), is in excellent 
agreement with the theoretical prediction ^ that includes thermal fluctua- 
tions. As S increases a small difference is observed. It is therefore likely that, 
as fluctuations become dominant, ^ is less accurate. It is therefore justified 
to employ Dynes ansatz to fit the experimental dl /dV in the nanoscale region 
provided that fluctuations are not dominant 5 <C max(Z\o, T c ). For larger fluc- 
tuations the generalized expressions ^ , ([6| provides a suitable generalization 
of Dynes expression in the nanoscale where now the fitting parameters are A 
and but I5 7^ Zo since _Td does not include the effect of thermal fluctu- 
ations. Therefore this generalization of the Dynes ansatz allows to single out 
the effect of fluctuations from other forms of decoherence. 



4 Quantifying fluctuations in STM experiments 

The analysis of experimental STM data by the SPA formalism (|5| , ^ devel- 
oped above has clear advantages over Q even for sizes in which Dynes ansatz 
describes well the energy gap. We focus on the study of r(T) which includes 
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Fig. 1 (a) dl/dV Q as a function of eV for S = 0.1A , T = 0.1Z\ , E D = 9A , A = 0.34. 
Red lines stand for the theoretical dl/dV from and dots correspond to Dynes fitting d5J. 
We assume p(e) fb 1/(5, (b) the order parameter for the same parameters as in (a). Red lines 
zl(T) roll correspond to the SPA prediction and black dots Au(T) stand for Dynes fitting, 
( c ) AhCw f° r t ne same parameters as in (b) (see text for more details). 



effects leading to a finite quasiparticle life-time and thermal fluctuations as 
well. We note that r = I/} in |2| englobes two contributions: l£> = r _1 + F t h 
where the former accounts for the quasiparticle life time and the latter is re- 
lated to thermal fluctuations arising due to finite size effects. It is therefore 
expected that r t h has a peak around the mean-field critical temperature with 
a typical width that describes the interval of temperatures where fluctuations 
are important. A important drawback of Dynes approach is that it is not pos- 
sible to disentangle these two contributions. The situation is different if the 
SPA formalism (|5|,([6]) is employed. Here thermal fluctuations are included 
so that 2~s = 1/r only accounts for the finite quasiparticle lifetime. This is 
an important difference. One can quantitatively estimate the role of thermal 
fluctuations by the quantity r t h = //) — Is which can be obtained by fitting 
the experimental data by both Dynes ansatz Q and the SPA expression (J5J. 
This quantity has a maximum (see Fig.JTJ:) near the critical temperature where 
thermal fluctuations are more important. Therefore it can be used to estimate: 
a) the importance of thermal fluctuations with respect to other sources of de- 
coherence, b) the (would-be) critical temperature and c) the region around 
T c where fluctuations are relevant. Similar findings are applicable to T <C T c 
provided that quantum fluctuations are also included in the theoretical formal- 
ism. In summary, the use of a theoretical framework that takes into account 
deviations from mean-field to analyse experimental data expands substantially 
the information that can be obtained from STM experiments. 



References 

1. R. C. Dynes, et al., Phys Rev. Lett. 53, 2437 (1984). 

2. G. Binnig, H. Rohrer, IBM Journal of Research and Development 30: 4 (1986). 

3. J. von Delft and D.C. Ralph, Phys. Rep., 345, 61 (2001). 

4. S. Bose, et al., Na ture Mat. 9, 550 ( 2010). 

5. I. Brihuega, et al., arXiv:0904.0354, Phys. Rev. B. in press. 

6. R. Denton, B. Mulschlegel, and D.J. Scalapino, Phys. Rev. B 7, 3589 (1973). 



